TGAL-92-15 


AD-A260  033 


SOURCE  MULTIPLICITY  EXAMINED  WITH 
MINIMUM  ENTROPY  DECONVOLUTION 


I.  H.  Henson  and  R.  K.  Cessaro 


Teledyne  Geotech  Alexandria  Laboratories 
314  Montgomery  Street 
Alexandria,  Virginia  22314-1581 


DTfC 

elects: 


NOVI  2  1992 


C 


OCTOBER  1992 


SEMI-ANNUAL  REPORT; 
ARPA  ORDER  NO.: 
PROJECT  TITLE: 
CONTRACT  NO.; 


No.  2  (4  April  1992  -22  October  1992) 

6731 

Multichannel  Minimum  Entropy  Deconvolution 
F29601-91-C-DB02 


Approved  for  Public  Release;  Distribution  Unlimited 


Prepared  for: 

PHILLIPS  LABORATORY 
KIRTLAND  AFB,  NM  87117-5320 


Monitored  by: 

DEFENSE  ADVANCED  RESEARCH  PROJECTS  AGENCY 
NUCLEAR  MONITORING  RESEARCH  OFFICE 
3701  NORTH  FAIRFAX  DRIVE 
ARLINGTON,  VA  22203-1714 


The  views  and  conclusions  contained  in  this  report  are  those  of  the  authors  and  should  not  be  interpreted  as 
representing  the  official  policies,  either  expressed  or  implied,  of  the  Defense  Advanced  Research  Projects 
Agency  or  the  U.S.  Government. 


9 


92-29286 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OM8  No.  0704-0188 


Public  reporting  burden  for  this  colleaion  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  ensting  data  sources, 
gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this 
collection  of  information,  including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  information  Operations  and  Reports,  12  is  Jefferson 
Oavis  Highway.  Suite  1204.  Arlington.  VA  22202-4302.  and  to  the  Office  of  Management  arsd  Budget.  Paperwork  Reduction  Project  (0704  0 108).  Washington,  DC  20S03 


1.  AGENCY  USE  ONLY  (Leave  blank) 


4.  TITLE  AND  SUBTITLE 


2.  REPORT  DATE 

22  October  1992 


3.  REPORT  TYPE  AND  OATES  COVERED 

Technical  Report,  4  April  1992  -  22  Oct  1992 


S.  FUNDING  NUMBERS 


Source  Multiplicity  Examined  with  Minimum  Entropy  Deconvolution 


6.  AUTHOR(S) 


Contract  F2960 1  -9 1 -C-DB02 


I.  H.  Henson  and  R.  K.  Cessaro 


7.  PERFORMING  ORGANIZATION  NAME(S)  ANO  AODRESS(ES} 

Teledyne  Geotech  Alexandria  Laboratory 
314  Montgomery  Street 
Alexandria,  VA  22314-1581 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


TGAL-92-15 


9.  SPONSORING /MONITORING  AGENCY  NAME(S)  AND  AOORESS(ES) 


Phillips  Laboratory  (PL/PKRC) 
Kirtland  AFB,  NM  87117-5320 


10.  SPONSORING /MONITORING 
AGENCY  REPORT  NUMBER 


DARPA-NMRO 

3701  N.  Fairfax  Drive  #717 

Arlington,  VA  22203-1714 


12a.  DISTRIBUTION  /  AVAILABILITY  STATEMENT 


12b.  DISTRIBUTION  CODE 


Approved  for  Public  Release;  Distribution  Unlimited 


13.  ABSTRACT  (Maximum  200  words) 


This  report  contains  the  preliminary  results  of  a  study  to  determine  the  usefulness  of 
minimum  entropy  deconvolution  in  detecting  seismic  source  multiplicity  and  its  potential  for 
discriminating  ripple-fired  explosions  from  other  seismic  events.  Specific  examples  of  its  appli¬ 
cation  to  quarry  blast  data  and  synthetic  seismograms  are  discussed. 


14.  SUBJECT  TERMS 


15  NUMBER  OF  PAGES 


source  multiplicity,  deconvolution 


16.  PRICE  CODE 


17.  SECURITY  CLASSIFICATION 
OF  REPORT 

Unclassified 


NSN  7540-01-280-5500 


18.  SECURITY  CLASSIFICATION 
OF  THIS  PAGE 

Unclassified 


19.  SECURITY  CLASSIFICATION 
OF  ABSTRACT 

Unclassified 


20.  LIMITATION  OF  ABSTRACT 


Standard  Form  298  (Rev  2-89) 

by  ANSI  Std  7)9- 10 

298  n; 


Teledyne  Geotecli 


F29601-91-C-DH02 


4  April  1992  -  22  October  1992 


1.  OBJECTIVES 

The  objective  of  this  project  in  the  first  year  is  to  develop  a  technique,  based  on 
minimum  entropy  deconvolution  (MED),  useful  for  discriminating  ripple-fired  explosions 
from  other  seismic  events. 

2.  PROGRESS 

During  the  last  quarter  we  experimented  further  with  data  from  a  known  quarry  blast 
and  synthetics.  Figure  1  shows  a  bore  hole  vertical  component  record  from  LNO,  for  a 
known  quarry  blast  at  a  distance  of  approximately  40  km.  After  some  experimentation 
with  the  filter  length,  damping  coefficient  and  data  window  position  and  length,  some  new 
results  have  been  obtained  for  this  data.  The  record  from  LNO  is  shown  as  the  top  trace  in 
Figure  1.  The  sampling  rate  was  60  sps.  The  middle  trace  shows  the  first  4  seconds  of  the 
signal  that  was  input  to  the  MED  algorithm.  One  4-second  MED  filter  was  generated 
which  deconvolved  this  signal  into  the  trace  shown  at  the  bottom  of  Figure  1.  It  is  tempt¬ 
ing  to  interpret  the  deconvolved  trace  as  4  phase  arrivals,  the  first  three  consisting  of  two 
closely  spaced  spikes  due  to  source  effects.  The  time  separation  between  the  spikes  in 
each  doublet  is  50  msec.  The  quarry  shot  layout  consisted  of  two  rows  of  shots,  27  shots 
total.  The  delay  time  between  individual  shots  was  approximately  9  msec.  The  delay  time 
from  the  start  of  the  first  row  to  the  start  of  the  second  row  was  42  ms.  Could  the  double 
nature  of  the  deconvolution  be  due  to  this  delay  time? 

Because  of  the  nonlinear  behavior  of  the  algorithm,  it  is  very  difficult  to  have 
confidence  in  such  an  interpretation  of  the  deconvolution.  Different  deconvolutions  can  be 
obtained  for  the  same  data  window  by  simply  varying  the  filter  length  or  damping 
coefficient.  Figure  2  shows  4  deconvolutions  for  the  data  window  of  Figure  1.  Each  trace 
in  Figure  2  was  generated  by  simply  varying  the  filter  length  or  the  damping  coefficient. 

We  have  tested  the  algorithm  on  synthetic  data  (60  sps)  generated  by  summing  time 
shifted  copies  of  real  data,  with  time  delays  equal  to  the  quarry  blast  delays  with  a  10  per¬ 
cent  random  perturbation.  Synthetics  were  also  made  with  simulated  time  delays  equal  to  a 
multiple  of  the  actual  time  delays.  Using  the  MED  algorithm  we  were  only  able  to  decon¬ 
volve  some  of  the  simulated  shot  times  when  the  delay  times  were  increased  by  a  factor  of 
5,  from  approximately  10  to  50  msec.  In  this  ca.se,  shown  in  Figure  3,  the  simulated  shot 
delay  time  was  approximately  equal  to  3  sample  intervals. 
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Figure  1.  top;  LNO/sz  recording  for  a  quarry-blast  at  a  distance  of  40  km.  middle;  the 
first  4  seconds  of  the  signal  used  as  input  to  the  MED  algorithm,  bottom;  one  of  the 
deconvolution  results  obtained  for  this  data  window. 
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3.  FUTURE  PLANS 

The  work  to  date  suggests  that  a  sample  rate,  higher  than  what  is  readily  available, 
will  be  required  in  order  to  recover  the  original  ripple-fired  quarry  blast  sequence  with  this 
method.  As  an  aid  in  investigating  the  relationship  between  ripple-fire  shot  spacing,  sam¬ 
ple  rate,  receiver  distance  and  recovery  efectiveness,  simple  synthetic  seismograms  with 
reasonable  shot  spacing  and  receiver  geometry  will  be  generated.  Once  the  limitations  of 
this  method  are  more  fully  explored  suitable  data  will  be  sought  and  examined. 


Deconvolutions  of  LNO/sz1  1 992042 


Figure  2.  Four  deconvolutions  obtained  for  the  data  window  shown  in  the  Figure  1.  Each 
trace  resulted  from  different  values  for  the  filter  length  and/or  damping  coefficient. 
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Figure  3.  Top:  Synthetic  generated  by  summing  data  27  copies  of  a  60  sps  LNO  record 
with  a  50  msec  time  delay.  The  synthetic  has  been  bandpassed  to  15-30  Hz.  Middle:  A 
deconvolution  of  the  synthetic  using  a  9  second  filter  Bottom:  The  deconvolution  shown 
at  an  enlarged  scale.  Approximately  9  of  the  27  shots  have  been  deconvolved. 
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8-CUBED.  LA  JOLLA 

MCLAUOFLN 

1 

M  WSTBL  OHCUTT,  GIVEN 

2 

STLOUtSU 

MrrCHaL  HERFMAW 

1 

SCLBED.  RE3TON 

MURfW 

z 

RADIX 

PULU 

1 

NORWAY.  NTNf' 

nNGOAL 

2 

TELEDYNE.  ALEXANDRIA,  VA 

RIVERS 

7 

TASC 

SAILOR 

1 

MAY-21-1991  08=25  FROM  DARPA 


TO 


915058469607 


COfiL  IIAUNG  UST4ai 


OtKIAMZATION 

MAIC 

Na  COPIES 

U8C 

SAUMIfi.  AKI 

1 

IRIS 

SMP9CN 

a 

UWASHMCnON 

SMrm 

1 

aWSOONSM 

THJRBmbCVe) 

1 

Mrr 

TGKSQZCANTY 

1 

U.AZ 

WAllACe 

1 

MOC.NEWWGTOKVA 

WORTMM4 

1 

US  OOVERNMENT  A08«CSS 

ACOA 

LOEPMAN 

1 

AFosnto* 

JBWPBW2X> 

1 

AfTAC.  C8S.  fWeSLVN.  VA 

BLANDfCRD 

1 

AFTACyCA 

sr»«%> 

1 

AFTAC/rr 

pajorre 

1 

CIA/AC(S 

KATISPOCEY 

1 

ClAOSWR 

TUWBUL 

1 

OARPA 

ALEWNERYALUKERR 

7 

QARPA/RMO 

LIBRARIAN 

1 

OIA 

GLCAH? 

1 

DNA/8PSS 

SK3RE 

1 

coe 

KDCNTZ 

1 

one 

KPOCTR 

2 

GIA.WH 

CIPAR 

1 

GIAWH 

LEWMdWCZ 

1 

GWXO 

X> 

1 

LlNl 

a 

OSD 

CORE 

1 

SAMCHA 

CRAB. 

1 

LB36 

LEnTI 

1 

U3S3 

UAS9E 

1 

wumBso 

RETKE 

1 

TOTAL  NUUBER  OF  REPORTS 

S8 

II 


